Supersolids in confined fermions on one-dimensional optical lattices 
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Using quantum Monte Carlo simulations, we show that density-density and pairing correlation 
functions of the one- dimensional attractive fermionic Hubbard model in a harmonic confinement 
potential are characterized by the anomalous dimension K p of a corresponding periodic system, 
and hence display quantum critical behavior. The corresponding fluctuations render the SU(2) 
symmetry breaking by the confining potential irrelevant, leading to structure form factors for both 
correlation functions that scale with the same exponent upon increasing the system size, thus giving 
rise to a (quasi)supersolid. 
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The recent experimental observation^ of nonclassical 
rotational inertia in solid 4 He revived the interest in the 
possible coexistence of superfluidity and crystalline order, 
a phase named a supersolid (SS), proposed long ago,— in 
spite of arguments against it 4 Such arguments seem to be 
validated by quantum Monte Carlo (QMC) simulations^ 
and theoretical work^ showing that a defect-free crys- 
tal of 4 He does not show off-diagonal long-range order 
(ODLRO). Further numerical simulations suggest that 
ODLRO may arise from domain walls^ or in a metastable 
state,— such that the existence of a pure supersolid phase 
in 4 He is still under debated 

Alternative paths toward the realization of a SS were 
recently suggested for quantum gases on optical lat- 
tices, e.g., for bosons with dipolar interactions^ Bosc- 
Fermi mixtures,— Bose-Bose mixtures,— frustrated lat- 
tice geometries^ or extended interactions 1 ^ by loading 
the atoms in higher bands Most of the proposals, how- 
ever, are based on mean-field approximations that tend 
to be unstable with respect to phase separation, as shown 
by QMC simulations^ Also the presence of a confining 
potential is in general neglected, discarding the coexis- 
tence of phase o 15 i 16 that may obscure the experimental 
determination of a SSJ^ 

Another system where a SS phase is well known to 
exists is the fermionic Hubbard model with an attrac- 
tive contact interaction [the first two terms in Eq. ([1]) 
below]. In contrast to the mentioned bosonic systems, 
where vacancies or defects lead to a SS, here diagonal 
and ODLRO result from an SU(2) symmetry, relating 
density and pairing amplitudes at density n = 1 J£ While 
the experimental realization of such a model is at least 
debatable in solid state systems, fermions with attractive 
contact interaction were recently achieved with quantum 
gases by means of Feshbach resonances, leading to pairing 
and superfluidity. 19 ' 20 i 21 Even the observation of super- 
fluidity with fermions in an optical lattice for a dense 
system (n ~ 1) was announced recently^ bringing very 
close the realization of an SU(2)-symmetric SS. 

Still, the presence of a confining potential could inhibit 
the formation of the SS, since it explicitly breaks the 
above-mentioned SU(2) symmetry, as discussed in detail 



below. We show, however, based on QMC simulations of 
the one-dimensional (ID) attractive fermionic Hubbard 
model in a harmonic potential, (i) that density-density 
(DD) and pairing correlation functions are characterized 
by the anomalous dimensional K p of the corresponding 
periodic model, and (ii) that in spite of the confinement, 
conditions can be reached where the structure form fac- 
tors of both correlation functions diverge with the same 
power as the system tends toward the thermodynamic 
limit. Hence, a proper treatment of long-ranged quan- 
tum fluctuations shows that the SU(2) symmetry break- 
ing becomes irrelevant, thus giving rise to a SS state. 
In the following, we consider the Hamiltonian 
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where cj CT and c J(T are creation and annihilation opera- 
tors, respectively, for fermions on site j (position Xj = ja, 
where a is the lattice constant) with spin a =T,4- The 
local density is nj = rijf + n,ji, where n J(T = c^c^. The 
contact interaction is attractive (U > 0) and the last term 
models the potential of the confining trap with TV (even) 
sites. The QMC simulations were performed using a pro- 
jector algorith m 24 i 25 ' 26 (T = 0) with a projector param- 
eter 6 ~ 30 /t that suffices to reach well-converged values 
of all observablcs discussed here. For time slices below 
At = 0.1/i, the corresponding systematic error has no 
effect on the results. We use in general At = 0.08/t. 
It is well knowni£ that a particle-hole transformation 
= (— iy Cj-f, dt^ = ct, applied to the Hamiltonian ((T|) 
leads to a repulsive Hubbard model in the presence of a 
spatially varying magnetic field hj = 2V(xj — Na/2) 2 . 
If (n iT ) = (n jU ), then n d } = (dtd^) + (d^d^) = 1, i.e., 
the new system is half filled. If, furthermore, hj = 0, i.e., 
for V = 0, a Mott insulator with dominating antiferro- 
magnetic correlations results for the repulsive system^ 
In terms of the operators of Eq. (1), this implies that 
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FIG. 1: Density (m) vs position Xi in units of the lattice 
constant a for Nf = 74 fermions in a harmonic confining 
potential Va 2 /t = 1.826 x 1CT 4 , with U/t = 4. The inset 
shows a detailed view of the density oscillations. 

for V = and n = 1 the DD correlation function 
Njg = (rijni) — (nj){ni) and the pairing correlation func- 
tion Pj£ = (A-Aj), with Aj = Cj^Cji, obey 

N 3l = 2(-l)^P ji , (2) 

and will be the dominating correlation functions, with 
the same power-law decay. This identifies the state as a 
SS. However, the presence of a confining potential breaks 
this SU(2) symmetry, so that for V^O the SS is expected 
to be destroyed. In the following we show that, due to 
quantum critical fluctuations, even for V ^ the SS is 
nevertheless recovered. 

A first important question is whether, in the presence 
of a confining potential, which produces a spatial de- 
pendence of the density, the correlation functions will 
each exhibit a decay with a single power. In fact, this 
is no longer the case if one relies on the local density 
approximation (LDA).— However, we find that the LDA 
is not appropriate to describe the physics of the model 
in Eq. (1). In Fig. Q] we show the density profile for 
Nf = 74 fermions, U/t = 4, and a characteristic density 
p = 1. As shown before,— the characteristic density de- 
fined as p = N ja{V It) 1 / 2 relates systems with different 
sizes, number of particles, and confining harmonic po- 
tentials in the same way as the particle density does for 
periodic systems of different sizes. Figure Q] shows that 
the presence of the confining potential produces density 
oscillations on a short distance scale (~ a), which are 
absent in periodic systems. These oscillations preclude 
the use of the LDA as a sensible approximation, since it 
would require a slowly varying density. 

We study next the DD and pairing correlations using 
QMC simulations. Since in the presence of quasi-long- 
range order a divergence is expected for a particular wave 
vector k in a periodic system, we examine the eigenvalue 
equation for DD correlations, 

^iV^=A^, (3) 
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FIG. 2: (Color online) Level plot of the moduli of the Fourier 
transforms of 0^ (a) and <fij (b) vs fi [see Eqs. © and Q]. 

which reduces to a Fourier transformation for a periodic 
system, i.e., the eigenvectors <pj are plane waves, each 
mode characterized by a particular k vector. Interest- 
ingly, we find that even inside a harmonic trap most of 
the modes <j>^ can still be assigned each to a characteristic 
wave vector. This can be seen by considering the moduli 
of the Fourier transform [~ y\ exp (ikxj)] of each . 
As shown in Fig. [Ha), most of the modes /i, with 
in ascending order, have a dominant component each of 
which can be related to a wave vector k (actually, both 
±fc due to reflection symmetry). We checked that re- 
gions where a clear identification of a dominant k vector 
is not possible at all correspond to the lowest eigenvalues 
that arise from the vanishing density at the edges of the 
fermionic cloud in Fig. [T] Results of a similar analysis 
for the eigenvectors Lp^ of the pairing correlations Pjg, 

£P^ = P<V?, (4) 

e 

are shown in Fig. HJb) . We find that the highest eigen- 
values of the pairing correlation function have the largest 
weight around k = 0, i.e., pairs are formed by fermions 
with opposite momenta. This unbiased result helps to 
understand why even in trapped systems, where trans- 
lation invariance is broken and momentum is not a 
"good" quantum number, pairing correlations have been 
observed predominantly between atoms with opposite 
momenta^ 
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FIG. 3: (Color online) JV(fc) for the confined (•) case with the 
same parameters as in Fig. [T] Superimposed are the results 
for a periodic system (□) with 2/cf given by the black vertical 
line. The resulting density is n = 0.71. Inset: position of the 
2&F peak as a function of U in the harmonic trap for p = 1. 
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FIG. 4: (Color online) Anomalous dimension K p extracted 
from the slope of N(k) for the confined system (x) vs U for 
p = 1.0 and (n) = 0.71. Superimposed are the results for a 
periodic system obtained from Bethe ansatz (©) for (n) = 0.7 
(Ref.[H). 



Since each mode p can be assigned rather clearly to a 
wave vector k, we can map each eigenvalue to N(k) 
as in a periodic system. Figure [3] shows the result of 
such an assignment (•), according to the salient features 
shown in Fig. [Ha). We superimposed in Fig. [3] the re- 
sults for a periodic system (□), for which the density 
was fixed by identifying the wave vector corresponding 
to the maximum in N(k) of the confined system with 
2k p of the periodic system (black vertical line), where 
kp is the Fermi momentum, establishing in this way a 
link between a trapped system with a characteristic den- 
sity p and a periodic system of density n = 2kpa/-n. Such 
an identification is based on the observation that, in a ID 
fermionic system with attractive interactions, 2kp oscil- 
lations in the charge channel are dominating. The inset 
shows the position of the peak as a function of U, which 
for the range 2 < U/t < 8 does not change appreciably. 
Hence, the identification is possible for a wide range of 
parameters and does not depend on any particular fine 
timing. 

While the identification of N(k) in the periodic and 
the confined systems deteriorates at wave vectors beyond 
2kp (with a number of points at zero or close to it due to 
the regions where the density vanishes), it is remarkably 
good close to k = 0. For a periodic system the DD cor- 
relation function obeys N(k) — > K p \k\a/ir for k — ► 0M- 
Therefore, the anomalous dimension K p that determines 
the power-law decay of the 2kp oscillations that dom- 
inate at large distances can be directly obtained from 
the slope of N(k) for k — ► 0. Figure [3] shows that the 
mapping proposed here leads to the identification of K p 
also for the confined system. In fact, such an identifica- 
tion is possible in general, as shown in Fig. 31 where a 
comparison of the values of K p obtained for a confined 
system with p — 1 and those obtained by Bethe ansatz 
for (n) = 0.7 is shown^ 

We are now in a position to determine the conditions 
to attain a SS. This is achieved in the periodic case for 



n = 1, i.e., 2kpa = n. We thus increase p until the 
maximum of N(k) in the confined system appears at 
k = 7r I a. This is shown in Fig. O where the mapping 
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□ periodic, U/t= 4, n=1.0 
• harmonic trap, U/t= 4, p=2.34 
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FIG. 5: (Color online) (a) N(k) for the confined case for 
p = 2.34 (•); other parameters as in Fig.fT] Superimposed are 
the results for a periodic system (□) with density n = 1.0. (b) 
P(k) for the confined (•) and periodic systems (□) with the 
same parameters as in (a). Inset: largest eigenvalues for the 
DD [N(k = 71")] and pairing [P(k = 0)] vs number of fermions 
Nf in a double-logarithmic scale. 

described above was made for both the DD [N(k)j and 
pairing [P(fc)] eigenvalues. Both quantities of the con- 
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fined system follow closely the behavior in the periodic 
system. As in the case of p = 1, the power-law behavior 
of the 2kp oscillations in N(k) for the trapped system 
is determined by the same K p as in the periodic case. 
Moreover, in the periodic case, due to the above men- 
tioned SU(2) symmetry, the relationship Eq. ((2]) holds 
between the correlation functions. This means in a peri- 
odic system that P{k) = N(ir/a — k)/2, a relation that 
is obeyed in Fig. [S] (□). Remarkably, and to a very good 
approximation, this symmetry is seen to be present also 
in the confined system (•). Hence, in spite of the ex- 
plicit breaking of the SU(2) symmetry by the confining 
potential, this symmetry is recovered in the correlations, 
which display quasi-long-range order. This fact is rem- 
iniscent of the universality found for hard-core bosons 
confined on a ID lattice, where the power-law decay of 
the one-particle density matrix remains unaffected by the 
presence of a confining potential^ In fact, in the limit 
U/t — ► oo of Eq. (1), the resulting on-site pairs reduce 
to hard-core bosons. Here, we find that, due to the long- 
ranged fluctuations, once the appropriate characteristic 
density is found, the SU(2) symmetry is restored with 
high accuracy for the considered values of U. 

The inset of Fig. [5] shows the scaling behavior of the 
largest eigenvalues N(k = if) and P(k = 0) for the har- 
monically confined system as a function of the number 
of particles for p = 2.34 and U/t = 4. Both indeed di- 
verge with the same power law, and, as discussed above, 



differ in magnitude within statistical errors by a factor 
of 2, as expected when the SU(2) symmetry is present. 
Hence, we show explicitly that quasi-long-range order is 
present in both the charge density and pairing channel 
with a unique exponent, giving rise to a SS state even in 
the trapped system. 

In summary, on the basis of QMC simulations (i) 
we have shown that quasi-long-range order in a con- 
fined system of attractive fermions is characterized by 
the single exponent K p of the related periodic system, 
and given a prescription to relate both systems; (ii) we 
have also shown that this identification leads to systems 
where diagonal and off-diagonal quasi-long-range order 
arise with the same exponent, such that they correspond 
to a supcrsolid in one dimension. These results clearly 
show that a SS state is attainable experimentally with 
fermionic quantum gases with two degenerate hypcrfine 
states confined in one-dimensional tubes with an axial 
optical lattice, a configuration that was already realized 
with fermions with repulsive interactions^ 

Note Added: After this work was completed and sub- 
mitted, we became aware of a related Letter by G. Xian- 
long et al£2- 
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